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Abstract 

We present bound state spectra of the 3D rational potential, V{r) = + Ar^/(1 + gr"^), g > 0, 
by means of the generalized pseudospectral method. All the thirty states corresponding to n=0-9 
are considered for the first time for a broad range of coupling parameters. These results surpass 
the accuracy of all other existing calculations published so far except the finite-difference method, 
which yields similar accuracy as ours. Variation of energies and radial distribution functions 
is followed with respect to the interaction parameters. Special emphasis has been laid on higher 
excitations and negative values of the interaction, where relatively less work has been reported. The 
energy sequence is found to be different for positive and negative interaction; numerically following 
a mirror-image relationship usually, if not always. Additionally, twenty energy splittings arising 
from certain levels belonging to n=0-9 are systematically studied as functions of the potential 
parameters. Several new states (including the higher ones) are presented. 
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I. INTRODUCTION 



Rational (also termed as non-polynomial oscillator, NPO) potentials are one of the most 
widely studied important model systems in quantum mechanics. The objective of this article 
is to investigate such potentials of the form, 

Ar2 



V{r) 



^2 _^ 



g>0 



(1) 



1 + gr'^ 

The Schrodinger equation (SE) can not be solved exactly in general and an impressive 
amount of theoretical works exist employing numerous formalisms as evident from a vast 
amount of literature [ll-j42|. It has significance in nonlinear Lagrangian methods in quantum 
field theory, laser physics, nonlinear optics, particle physics, etc. jl-^. Obviously for A = 
orA = (7=Oor(7<<Aor large g, its solution behaves as the harmonic oscillator. The 
potential in ID has attracted enormous attention from a large number of authors for nearly 
three decades employing several powerful tools such as variational method, perturbation 
theory (PT), semi- numerical or purely numerical technique. Exact analytical solutions of 
ground- and some excited-state eigenvalues and eigenfunctions in ID are possible provided 
g and A are related in certain specific ways j^y. 
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Some of the most important 



works are chronologically: Ritz variational principle using Hermite polynomial basis 



12|, 



variational and PT calculations 13|], PT using [6,6] Pade approximation and hypervirial 



theorem [l^, a variety of finite-difference approaches o: 
PT Q, 

combined supersymmetry and WKB approach 



different flavour 



ISHlTJ, algebraic 



111 ], continued fraction method 22| . 



analytic continuation method [23], PT with mixed hypervirial and Hellmann-Feynmann 



;heory 



30|, Hill-determinant method |3l|, variational bounds via Rayleigh-Ritz theorem 



321 ]. a quadrature discretization technique [33[ and purely numerical calculation [35|, etc. 
Some of these deal with only ID and others for both ID plus 3D and/or N dimensions. 

In parallel to the ID case, 3D NPO has also been treated in significant detail in the 
past and recent years, and the interest continues to grow. Exact solutions were reported for 



some special cases of g and A using SUSY 



generalized harmonic oscillator form 



3, 



etc. However since these are of limited applicability, several 



Some of the promising ones are: shifted 1/N expansion 38 



methodologies were attempted. 



39| , bounds from eigenvalue mo- 



ments method 40j, hypervirial calculation 411 , bounds through the envelope method 42]. 
However as clear, 3D NPO remains relatively less explored as compared to its ID counter- 
part. Amongst these, satisfactory results for arbitrary combination of potential parameters 



of a general state can be obtained only by a very few. Thus while for ID case, more than 
decimal-place accuracy eigenvalues were reported in a number of papers (for example 
33| . to our knowledge, only a few of the mentioned works 4CH42j are capable 
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23 



30 



of delivering such results for 3D NPO satisfying all these criteria. Moreover two of these 
works in |40| and 42| deal mostly with estimate of the bounds, not the direct values. Thus 



there is a scarcity of good-quality results for 3D NPO as compared to that of ID. It may 



be mentioned here that with the rare exception of 



practically all works on 3D NPO 



have focused on (+)ve A, even though it has been known for long time that A<0 also offers 



31 



well-behaved solutions provided 5'>0. Negative A has been examined in ID in detail 
however. 

In recent years, generalized pseudospectral (GPS) method has witnessed great success 
and been proved to be a powerful, viable alternative for accurate and efficient treatment 
of a variety of potentials in quantum mechanics, such as the spiked harmonic oscillator, 
the power-law, the logarithmic, the Hulthen, Yukawa, etc., as well as both static, dynamic 
properties of many-electron systems including multiply excited high-lying Rydberg states of 



atoms [43|-[50|]. Successful applications have also been made for the 3D spherical confinernent 
studies of the isotropic harmonic oscillator, the H atom and the Davidson oscillator [50]. A 
review of the GPS method in the context of atomic and molecular calculations is available in 



471]. For some of the problems such as those in refs. 



45 



46 



48 



49| . this has offered results 



either quite comparable to those obtained from the best available methods in the literature 
or has surpassed the existing calculations. Another attractive feature is that it was very 
successful to furnish equally high-quality results for both lower and higher states, as well as 
for weaker and stronger interaction, which was not possible for many of the commonly used 
methodologies (see, for example 49|). However, to our knowledge, the GPS method has not 
been applied in the context of NPO so far. Hence our first objective is to extend its scope 
and validity in the realm of this heavily studied important potential. As already hinted, we 
will treat both A > and A < with special attention to the higher-lying states. To this 
end, calculated energies are presented for a wide range of g, A. Then the changes in energy 
are monitored with respect to the parameters. Secondly we explore the possible patterns 
in energy orderings by systematically performing all the thirty states with varying Ur and 
/ quantum numbers lying within n=0-9 for both (+)ve and (— )ve A. Besides, it is known 
that the characteristic degeneracy of certain energy levels of the 3D harmonic oscillator is 
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removed in certain situations (in the limit of A— J-O or q << A, for example); first four of such 



splittings have been discussed at length for (+)ve A j22|, |38| corresponding to n=0-4. Here 
we extend such an investigation for all the twenty differences possible within the range of 
n = — 9 for both (+)ve and (— )ve A as functions of the coupling parameters. Comparison 
with literature data has been made wherever possible. Section II gives a brief account of 
the method of calculation. Discussion on the results are given in Section III while a few 
concluding remarks are made in Section IV. 

II. METHOD OF CALCULATION 

This section summarizes the essential details of GPS formalism to solve the nonrelativistic 
radial S 



IpnA'^) = ^n/ i^nli'^) (2) 



of a single-particle Hamiltonian. Relevant details have been presented elsewhere 
43|]-|50| and references therein); thus not repeated. Unless otherwise mentioned, atomic 
units are employed throughout this article. 

The desired SE can be written in following form, 

1 d2 £{£+!] 
h -^^ + v(r) 

2 dr2 ^ 2r2 ^ ^ ' 

where v{r) = V{r)/2. A 1/2 factor is introduced for sake of consistency with the literature. 

The key step is to approximate a function f{x) defined in the interval a; G [—1,1] by an 
N-th order polynomial /Ar(x) such that, 

f{x) = h{x) = j2f{x,)g,{x), (3) 

j=0 

ensuring that the approximation is exact at the collocation points Xj, i.e., fN{xj) = f{xj). In 
the Legendre pseudospectral method used here, xo=—l, XAr=l, while Xj{j = 1, . . . , A^ — 1) are 
obtained from the roots of first derivative of the Legendre polynomial Pn{x) with respect 
to X, i.e., P']y{xj) =0. The cardinal functions gj{x) given by, 

N{N + 1)Pn{x,) x-x, ' 
have the unique property, gj{xj>) =Sjij. Now (i) mapping the semi-infinite domain r G [0, oo] 
onto the finite domain xG[— 1,1] by the transformation r = r{x), (ii) using an algebraic 
nonlinear mapping, 

r = r{x) = L , (5) 

1 — X + a 

(iii) followed by a symmetrization procedure leads to the following equation 



TABLE I: Calculated lowest (rir = 0) eigenvalues E (times 2 in a.u.) of the 3D NPO for several q 
and A along with literature data for i = 0, 1,2. Asterisks denote exact analytical values [39|, I42l |. 

fm.nx — 150 a.u. 
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-133.000000000000 


-133.000000000000", 






















-133* 



42| 



"Variational calculation 
''Shifted 1/N expansion result 



AT r 



^l^j'j + v{r{xj)) + 6j>j Vm{r{xj)) 



A, 



EAj>, 



A^- 1, 



(6) 



where Aj = [r'{xj)]^^'^ ip{r(xj)) [Pjv(a;j)] ^; symmetrized second derivatives of the cardinal 

(2) 

function, Dj,j are given in the references [43|-|49|. 

By performing a series of test calculations, a consistent set of parameters were chosen 
(a=25, N=200) which produced "stable" converged results. These were used for all the 
calculations reported in this work. Usually Vmax was set at 150, but for higher excitations 
larger values up to 300 a.u. was used (see tables and Section III for details). The energies 
are given only up to the precision that maintained stability and are truncated rather than 
rounded-off. 

III. RESULTS AND DISCUSSION 



At first, we examine the convergence of our calculated energy eigenvalues of the 3D NPO. 
Table I compares a number of levels for particular values of g, A which offer exact analytical 
results (denoted by asterisks). These are available only for A<0 and given for the lowest states 
{nr=0) having /=0, 1, 2. Note that all energies in this table and throughout the article give 
E multiplied by a 2 factor, for easy comparison with literature. In all cases, 12-13 decimal- 
place accuracy is easily obtained. The variational calculations of 42| employ a Gol'dman 
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TABLE II: First two {rir = 0, 1) eigenvalues E (times 2 in a.u.) corresponding to / = — 3 of 3D 
NPO for select g and A along with the literature data, r^ax = 150 a.u. 













E 




9 


A 


This work 


Reference 





0.1 



0.1 



100 



100 



3.120081864016 
7.231009980656 
5.186373002931 
9.276541985488 

7.2439618404219 

11.317997742355 
9.294359110874 

13.355727291254 
3.507388348905 
7.648201241719 
5.651393306756 
9.713754138848 
7.734828038042 
11.76582837889 
9.787669778003 

13.804700633187 
3.983098339488 
7.984443523273 
5.993438790399 
9.993516159965 
7.996024670900 

11.996039234884 
9.997153638476 
13.99715862578 



3.1200" 
7.2312" 

5.1863730029314<E<5.1863730029316\ 5. 1863730029314^5. 186338'',5. 1864" 
9.276635'^,9.2766" 

7.2439618404138<E<7.2439618404260^ 7.2439618404189S7.243927'', 7.244" 
11.258" 

9.29435911086337<E<9.29435911088159', 9.2943591108746^9.2944" 

3.50738835'', 3.50738835%3.5122" 
7.64820124^ 7.64820124=, 7.6252" 

5.6503<E<5.6521*, 5.651393317250^5.652112'',5.6522" 
9. 705584 '',9.7056" 

7.734<E<7.736^ 7.734828042923^7.734778'^, 7.7348" 
11.7630" 

9.7875<E<9.7881', 9.787669778509^9.7876" 

3.983098^,3.9844" 
7.984444.'',7.9910" 

- <E<6.389\ 5.993438873366S5.993439A5.993565'^,5.9936" 
9.993516J",9.994694'', 9.9946" 

7.9947<E<8.037800*, 7.996024673021S7.996025A7.996048'', 7.9960" 
11.996039-',11.9964" 

9.9969<E<10.0113', 9.997153638602^9.997145^9.9972" 



"Shifted 1/N expansion result |38|| . "^Quadratic discretization rijisult jSSj. ^Continued fraction result 

oTh eigenvalue moment method [40|. 



''Lower and upper bounds £ 
"^Variational upper bounds 
"^Shifted 1/N expansion result 



42| 



and Krivchenkov Hamiltonian as a solvable model and there is excellent agreement of our 
results with theirs for all cases. Some of these were reported long times ago through shifted 



1/N expansion [39| and quoted appropriately. The present results clearly outperform them. 
As a second test, a few representative calculations were performed with A=0 for several 
low and high /. In all occasions, expectedly harmonic oscillator eigenvalues were recovered 
promptly to at least 13th place of decimal and hence omitted. This amply demonstrates the 
power and authenticity of GPS method in the present context. 



Now we are ready to present our central result in Table II. The calculated eigenvalues 
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g=100 


/■' /' /g=o-i 


9=10 


9=1/ 


'9=0.5 1 





FIG. 1: Variation of the ground-state energy (2E) of 3D NPO with respect to (a) g for fixed A, (b) 
A for fixed g. 

are compared with the best existing hterature data for three representative pairs of {g, 
A), viz., (0.1,0.1), (1,1) and (100,100), denoted respectively as (a)-(c), covering sufficiently 
broad ranges of interaction. For each of these sets, lowest two values (^^=0,1) belon ging 



to /=0, 1,2,3 are given. Bounds were estimated from the eigenvalue moment method [40| 
for the first states belonging to / >1 of all them. As seen, for smaller values of g and A in 
(a), the bounds are quite good; but as one passes to (c) for larger values, they deteriorate 
quite badly. Significantly improved upper bounds for same states ((b) and (c)) have been 
published recently For the highest state considered (?t,j.=1,/=3) of each pair (a)-(c), no 
results could be found for comparison. The shifted 1 /N expansion results are reported for all 
of these states except the highest ones. The first seven states of (c) have also been studied 



by a continued fraction method j22|] with modest accuracy. Clearly our results turn out to 
be the most accurate direct estimates for all of these cases. Now Fig. 1 displays the changes 
in ground-state energy (times 2) against g for seven fixed values of A (—100,-50,-10,1,100, 
200) in (a) and A for five fixed values of g (0.1,0.5,1,10,100) in (b). Similar plots hold 
for higher energy levels. For a fixed A, with increase in g, eigenvalues decrease steadily 
approaching those of the harmonic oscillator asymptotically for large g, as expected, and 
a larger A shows this behavior with greater magnitude. Correspondingly opposite trend is 
noticed for negative A. In (b), we see that the changes against A are more prominent for 
smaller g and as g increases, again eigenvalues approaching those of the harmonic oscillator. 
For further test, in Table III, we show some specimen results of the 3D NPO for some 
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TABLE III: High- lying eigenvalues E (times 2 in a.u.) of the 3D NPO potential for some selected 



g and A along 


with the literature data, rmax = 300 


a.u. 










E 








E 


i 9 


A 


This work 


Ref. [41] 


9 A 


This work 


Ref. [41] 


10 0.1 


200 


311.860880892760" 


311.860880892760 


0.2 1000 


699.105622574512* 


699.105622574512 






361.308788624591 


361.308788624592 




811.026062535266 


811.026062535263 






409.563564997653 


409.563564997653 




920.708010354176 


920.708010354177 






456.909323126874 


456.909323126874 




1028.176013499333 


1028.176013499331 






503.090316800161 






1135.454829499741 




20 0.5 


500 


914.366310994354"= 


914.36631099435 


0.1 1000 


1312.25167480938'' 


1312.251674809389 






990.80662152100 


990.80662152100 




1425.58600163865 


1425.586001638665 






1066.14830339382 


1066.14830339382 




1537.79045903952 


1537.790459039524 






1140.40019428154 


1140.40019428154 




1648.87112288517 


1648.871122885173 






1213.57118403531 






1758.83409624488 





"Lower 
^Lower 
'^Lower 
''Lower 



and upper bounds 
and upper bounds 
and upper bounds 
and upper bounds 



3n.8601371<E<311.8616 
699.10424<E<699.10909 
914.36540<E<914.36851 
1312. 25006<E<1312. 25333 




higher excited states. For this we choose two sets of {g,X)', for each of them first five states 
befonging to Z=10,20 are considered. In the literature, very few results are available for large 
A as this often creates problems and this motivates us to include higher A. The high- lying 
excited states are diffuse and extend over a large spatial region. Thus to incorporate these 
long-range contributions properly, larger values (300 a.u.) of Vmax needed in these cases. 
Lower and upper bounds for the lowest state of each set were reported [40| and current 
results are in consonance with those estimates. First four states were also calculated by 



finite-difference method j4l| and as seen, GPS results virtually coincide with those accurate 
values. In many cases, these two offer completely identical energies; for others they differ in 
the last decimal place. For nr=4, no reference could be found. 

Table IV collects eigenvalues for four sets of coupling parameters; two for (+)ve and 
)ve A each. This table serves three purposes: (i) give energies for A <0, for which no 



calculations have been made except those in [SS], (ii) give similar results for A > 0, for 
which, as already mentioned earlier, accurate estimates, especially for higher energy levels, 
are quite scarce, and (iii) make a comparative study to discern the energy orderings in 
(+)ve and (— )ve A regimes side by side. To accomplish this, we systematically calculated 
all the 30 energy levels of n ranging from to 9 for forty-five sets of {g,X), with g values 
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TABLE IV: Calculated eigenvalues E (times 2 in a.u.) of 3D NPO for four pairs of g and A in 



parentheses along with literature data for all the states corresponding to n=0-9. r^ax = 150 a.u. 





n 


(1000,0.1) 


(10,1000)"''' 




71 


(0.1,-1)" 


(1,-100) 


Is 





3.00009981081 


64.8250831107 


Is 





1.5624896124 


-79.086464632 


Ip 


1 


5.00009993345 


89.1234516907 


Ip 


1 


2.8620942889 


-78.861621366 


2s 


2 


7.00009981589 


94.8759670627 


Id 


2 


4.2817900651 


-78.421433572 


Id 




7.00009996001 


100.703995819 


2s 




4.4940422546 


-75.296619703 


2p 


3 


9.00009993351 


101.225823801 


1/ 


3 


5.7883629695 


-77.782769536 


1/ 




9.00009997141 


105.507579307 


2p 




6.0427339658 


-75.063364459 


3s 


4 


11.0000998195 


103.184601259 


Iff 


4 


7.3613625227 


-76.966703116 


2d 




11.0000999599 


105.945405717 


2d 




7.6312350478 


-74.608167024 


l9 




11.0000999777 


108.527833520 


3s 




7.7681684221 


-71.510476611 


3p 


5 


13.0000999335 


106.763936963 


Ih 


5 


8.9870011416 


-75.995560432 


2/ 




13.0000999713 


109.683699997 


2/ 




9.2593659985 


-73.950607078 


Ih 




13.0000999817 


111.060023667 


3p 




9.4398078008 


-71.267997180 


5p 


9 


21.0000999336 


115.948021680 


11 


9 


15.851969149 


-70.958113492 


4/ 




21.0000999712 


117.926235639 


2i 




16.091742759 


-69.765978162 


3h 




21.0000999816 


119.098277713 


3h 




16.309144680 


-68.250662394 


2j 




21.0000999864 


119.632225052 


4/ 




16.483741176 


-66.287773320 


U 




21.0000999893 


119.930276951 


5p 




16.592021371 


-63.686362266 



"First nine levels are: 64.825083, 89.123452, 94.875967, 100.703996, 101.225824, 105.507579, 103.184601, 



105.945406, 108.527834, ref. [2^. 
^First nine levels are: 64.740, 89.040, 93.360, 100.726, 101.108, 105.496, 99.440, 106.030, 108.528 
"^First nine levels are: 1.5604, 2.8614, 4.2814, 4.5098, 5.7882, 6.0512, 7.3612, 7.6364, 7.8164, ref. 



raf. 



0.1,1,10,100,1000 while the latter as -100,-10,-1,-10.1,0.1,1,10,100,1000. Only seventeen 
of these states (corresponding to n = — 5,9) are presented here, which is sufficient to 
support the conclusions drawn herein. Note that in this table, we follow spectroscopic 
notation, i.e., the levels are labeled as n^+l and / values. So 72^=3, 1=2 signifies Ad level 
and so on. Out of all these {g,\) pairs, four are presented for reasons to become clear in 



the following. It may be mentioned that in 



22 



energies up to n=4 leading to nine 



combinations of n^,/ {Is-lg) were given for (+)ve A, although in a different context (see 
later). For (+)ve A, we find that the ordering given in column 1 holds true for twenty-three 
sets out of twenty-five, excepting (1,100) and (10,1000). That is why (10,1000) is picked 
up. The first instance of such a violation is encountered for between the levels (1/, 3s) for 



A > 0, and (2s, 1/) for A < 0. The first nine states of the set (10,1000) reported in [22|, l38| 
interestingly show this feature qualitatively even though present values are far superior to 
those. Thereafter this ordering is not followed in many occasions as can be seen from the 
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table. The (1000,0.1) set, besides being a representative of the ((?, A) set obeying the most 
observed sequence, also numerically illustrates that as g increases, eigenvalues approach the 
harmonic oscillator values and in the limit of A— J-O, all the levels of a particular n tend to be 
degenerate. Similar observations also hold for (1000,-0.1), except that now the eigenvalues 
approach harmonic oscillator values from below. Out of twenty {g-,X) pairs with (— )ve A, 
the ordering in column 5 is found to be valid for seventeen sets while the sets (0.1—10), 
(0.1,-100), (1,-100) do not follow this. So (0.1,-1) is a representative of those seventeen 
sets and also in this case first nine levels have been reported earlier in the literature (given 
in footnote). The (1,-100) set does not follow the same trend as (0.1, —1) and is given in 
last column. We note that the "usual" or mostly observed ordering in these two cases of 
(+)ve and (— )ve A for a given n are mirror images of each other (or reversed). For example, 
when ?T,=9, the energy ordering is 5p< 4/ < 3/i< 2j < 1/ for (+)ve A, while for (— )ve A it is 
l/<2j<3/i<4/<5p 

Now Fig. 2 displays the changes in radial distribution for the first two states (1/ and 3s) 
in the vicinity of g and A values, where the energy ordering is broken. The bottom panel 
((a), (b)) corresponds to If state, while the top panel ((c), (d)) represents 3s state, (a) 
and (c) show the density changes for a fixed A (1000) while varying g] (c) and (d) show the 
same for a fixed g (10) with respect to variations in A. For both the states, we note that, 
for a fixed A, with an increase in g, the peak positions shift to the higher values of r, the 
peak height decreases and the distribution broadens. For a fixed g, with an increase in A, 
correspondingly opposite trend is observed in (b) and (d) for both the states. 

It is well-known that for a 3D harmonic oscillator ra^ and / quantum numbers satisfying 
n = 2nr + I are degenerate. Thus for n=4, (n^,/) pairs having values (2,0), (1,2) and (0,4) 
are degenerate, etc. For non-zero A, such degeneracies vanish and these are conveniently 
analyzed through the respective level spacings, AE = ; — E„,/ //. Quite detailed studies 
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We found similar 



were conducted for first four (n<A) such splittings of positive A 
qualitative features for those in 22| and thereby not repeated here. However, to our knowl- 
edge, no such attempts have been made as yet for A <0. Table V gives variations of all the 
first twelve such splittings possible between certain adjacent levels of n=2-l with respect to 
the interaction parameters. First two splittings Ei q — Eo,2 and Ei i — Eq^s belong to n=2,3 
respectively, while the last three Ei^s — Eqj, E2,3 — Ei^s, Es^i — E2,3 correspond to n=7. In 
columns 2-5, g is kept fixed at 0.1, while A varied from —0.1 to —100. Clearly all these split- 
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FIG. 2: Variation of the radial probability distribution functions, |ri?„^£p, for the 1/ (bottom 
panel) and 3s (top panel) states of the 3D NPO. For (a) and (c), A = 1000; for (b) and (d), g=10. 
See text for details. 



tings tend to increase as |A| increases; a trend as also observed for A>0 [22']. Furthermore, 
according to them the first four splittings, i.e., the reverse of ours (Eo,2 — Ei^o, Eo,3 — Ei^i, 
El 2— E2,o, Eo,4— El 2) increases continuously with A, eventually approaching a constant value 
in the limit of A — ?■ 00. However they considered A > 10^ and more elaborate calculations 
would be required to confirm whether such a trend will also occur for A<0. Similar trend is 
observed for other g as well as the remaining eight splittings (not given) occurring between 
certain permissible levels within n=% and 9. Finally we examine the changes in spacings 
with respect to g (from 0.1-1000), keeping A fixed at —100 (in columns 5-9). Thus the 
entries belonging to (0.1,-100) in column 5 are common to both the variations and hence 
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TABLE V: Variation of first twelve splittings (AE) of 3D NPO occurring between certain adjacent 
levels of n=2-7 with respect to A (left), and g (right). Parentheses enclose the respective {g,X) 



pairs. See text for details. 



No. 


(0.1,-0.1) 


(0.1,-1) 


(0.1,-10) 


(0.1,-100) 


(1,-100) 


(10,-100) 


(100,-100) 


(1000,-100) 


1 


0.014016 


0.212252 


3.008518 


3.727614 


3.124814 


0.528414 


0.011480 


0.000144 


2 


0.018938 


0.254371 


2.828010 


3.683210 


2.719405 


0.252782 


0.003626 


0.000038 


3 


0.021770 


0.269873 


2.657399 


3.638930 


2.358536 


0.142225 


0.001742 


0.000018 


4 


0.010083 


0.136933 


2.976420 


3.727042 


3.097690 


0.446619 


0.010609 


0.000140 


5 


0.023245 


0.272365 


2.497516 


3.594793 


2.044953 


0.089829 


0.001026 


0.000010 


6 


0.014199 


0.180442 


2.786802 


3.682481 


2.682610 


0.232326 


0.003561 


0.000038 


7 


0.023833 


0.268071 


2.348674 


3.550816 


1.776512 


0.061534 


0.000678 


0.000001 


8 


0.016888 


0.203685 


2.609170 


3.638047 


2.317747 


0.135666 


0.001730 


0.000018 


9 


0.007716 


0.098492 


2.940544 


3.726466 


3.068309 


0.384311 


0.009937 


0.000137 


10 


0.023845 


0.260169 


2.210768 


3.507016 


1.548542 


0.044682 


0.000482 


0.000005 


11 


0.018550 


0.214770 


2.444401 


3.593757 


2.004532 


0.087238 


0.001023 


0.000010 


12 


0.011165 


0.136694 


2.740799 


3.681748 


2.643349 


0.214832 


0.003500 


0.000038 



enclosed in box with border from both sides. Clearly they tend to disappear in the limit of 



large g, in agreement with a conclusion found in 22 



for positive A. Once again, this 



has been found to be numerically true for the rest eight splittings occurring between higher 
levels (not presented in the table) as well as for other values of A (—0.1,-1,-10). 

Before passing, it may be noted that, in the last decade, the power series solution method 
(originally proposed by 5l|, l52]) and some of its variants have been very successful in pro- 



ducing quite accurate eigenvalues, radial expectation values, etc., of a number of systems 
such as the confined 3D isotropic harmonic oscillator, 2D H atom confined in a circle with 



impenetrable walls, 3D H atom confined in an impenetrable sphere |53l-l55l|. Very recently, 
highly accurate eigenvalues have also been reported 55[] by means of a formal solution of 
the confluent hypergeometric function for the confined H atom. It may be interesting to 
employ such methods in the context of the NPOs. 



IV. CONCLUSION 



Quantum mechanics has nowadays widely spread applications almost everywhere in con- 
temporary science and technology. However, exactly solvable quantum systems are very 
scarce. Nearly exact numerical solutions of quantum-mechanical problems is of paramount 
interest in this vein, and also for the merits of the theory itself. This work offers nearly ex- 
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act solutions for a specific class of quantum potentials, the non-polynomial oscillator class, 
that has broad applications in various quantum and quantum-field models. Accurate bound 
states of the 3D NPO are calculated and considered in detail. While several high-quality 
results have been published for the ID NPO, the data for 3D NPO is quite limited. Reliable 
variational bounds have become available only very recently. Moreover, the A<0 region is 
explored very little. For both positive and negative A, as well as for far high-lying states, the 
present GPS calculations match very well with the best literature results. All levels in the 
range of n=0-9 (comprising thirty overall) are investigated for a wide range of the coupling 
parameters. This study surpasses all the existing results except the finite-difference data of 
ref. 4l|, which offers similar accuracy. 



Energy and radial density distribution variations are followed as functions of g and A. An 
attempt is made to identify the energy orderings for A <0 and A >0 separately. Although 
there seems to be a qualitative pattern, several deviations of oscillatory character to it 
are observed. Finally all the twenty energy differences between certain admissible levels for 
n=0-9 are examined for the first time as functions of the parameters. The method employed 
here is simple, easily applies to a wide variety of strong and weak interactions and performs 
very well for low and higher excitations. We hope that our results may serve as benchmarks 
for these potentials. 
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